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Abstract Thermal fluctuations of the color superconducting order parameter in dense 
quark matter are investigated in terms of the phenomenological Ginzburg - Lan- 
dau approach. Our estimates show that fluctuations of the di-quark gap may 
strongly affect some of thermodynamic quantities even far below and above the 
critical temperature. If the critical temperature Tc of the di-quark phase transi- 
tion were rather high (;S (50 -r- 70) MeV) one could expect a manifestation of 
fluctuations of the di-quark gap in the course of heavy ion collisions (above Tc). 
For Tc ~ 50 MeV color superconducting fluctuations may also affect an initial 
stage of the hybrid star evolution. 
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Introduction 

The quark-quark interaction in the color antitriplet channel is attractive driv- 
ing the pairing, cf. [1]. The problem has been re-investigated in a series of 
papers following Refs. [2, 3], see review [4] and Refs therein. The attraction 
comes from the one-gluon exchange, or from a nonperturbative 4-point inter- 
action motivated by instantons [5], or nonperturbative gluon propagators [6]. 
The zero-temperature pairing gap A was predicted to be ~ (20 200) MeV 
for the quark chemical potentials fig ~ (300 500) MeV. In the standard BCS 
theory, cf. [4], the critical temperature is estimated as Tc ~ 0.57A. 

One expects the di-quark condensate to dominate the physics at densities 
beyond the deconfinement/chiral restoration transition and below the critical 
temperature. Various phases are possible. E.g., the so called 2-color supercon- 
ductivity (2SC) phase allows for unpaired quarks of one color. There may also 
exist a color-flavor locked (CFL) phase [7] for not too large value of the strange 
quark mass m^, for 2A > ml/ fig, cf. [8], where the color superconductivity 
(CSC) is complete in the sense that the di-quark condensation produces a gap 
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for quarks of all three colors and flavors. The values of the gap are of the 
same order of magnitude for 2SC and CFL phases, whereas relations between 
critical temperature and the gap might be different, 0.57A for 2SC and 
Tc ~ 0.7A for CFL phase [9]. There are also another possibilities, e.g., of 
the pairing in spin-one channel, for which the pairing gap proves to be small 
A < 1 MeV, see [9]. 

The high-density phases of QCD at low temperatures may exist in interi- 
ors of most massive neutron stars (so called hybrid stars containing the quark 
core and the hadron shell) affecting the cooUng, rotation and magnetic field, 
cf. [10-12]. It is also possible to ask whether CSC is relevant for the terres- 
trial experiments? To produce a color superconducting matter in the laboratory 
one would need to cook a dense and not too hot baryon enriched matter. The 
nuclear matter prepared in heavy ion collisions at SIS, AGS, SPS or RHIC 
in all the cases has presumably not sufficiently high baryon density and, on 
the other hand, the matter prepared at SPS or RHIC is feasibly too hot in or- 
der to expect a manifestation of the CSC. The most relevant is probably the 
GSI "Compressed baryon matter" future heavy ion collision facility which may 
cook enough dense and not too hot state. Baryon densities up to ten normal 
nuclear matter density Wpo) at temperatures T < 170 MeV are expected 
to be reached at an initial colfision stage. It is supposed that the system is in 
the quark-gluon plasma state at such conditions. Then the system expands and 
cools down. In this process the temperature decreases up to T ~ 140 MeV 
at still rather high baryon density at an intermediate colfision stage. Although 
the temperatures at relevant densities are most likely larger than the critical 
temperature Tc of CSC, one may rise the question on a possible manifesta- 
tion of the precursor phenomena of the CSC phase transition, if Tc is rather 
high (Tc ^ (50 70) MeV). Recently, ref. [13] considered such a possibility 
within the Nambu-Jona-Lasinio model and demonstrated that the fluctuating 
pair field results in a prominent peak of the spectral function, which survives 
in the temperature interval |r — Td ^ (0.1 0.2)rc. Besides, it is interesting 
to investigate the role of the order parameter fluctuations for T < Tc, if Tc 
is ~ 50 MeV, that still may affect the neutrino radiation of the most massive 
hot neutron stars (if they indeed have quark cores) and the heat transport at an 
initial stage of their evolution. 

This paper is an extended version of the work [14]. We will study precursor 
phenomena of the CSC in the framework of the phenomenological Ginzburg - 
Landau approach. To be specific, we will consider condensates with total an- 
gular momentum J = that are antisymmetric in color and flavor. Such pair- 
ing states can occur in the weak coupling limit because one-gluon exchange is 
attractive in the color anti-triplet channel. Although for reasonable values of 
densities and temperatures conditions of the applicability of the weak coupUng 
Umit are hardly fulfilled, in order to get a feefing on possible relevance of flue- 
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tuation effects we will still use this limit. At the same time we will by hand 
vary the relation A(/iq), since the corresponding exponential dependence is 
most sensitive to corrections of the running QCD coupling constant compared 
to its perturbative value. We will restrict our discussion by the consideration 
of phases, which feasibly have large gaps, Uke 2SC and CFL, as the most in- 
teresting case for applications. 

Physics of pairing fluctuations 

Like there always exists a vapor under the water, there are excitations on the 
ground of any condensate. They appear due to quantum and thermal fluctua- 
tions. In classical systems and also at not too small temperatures in quantum 
systems, quantum fluctuations are suppressed compared to thermal fluctua- 
tions. Excitations are produced and dissolved with the time passage, although 
the mean number of them is fixed at given temperature. Pairing fluctuations 
are associated with formation and breaking of excitations of a particular type, 
Cooper pairs out of the condensate. Fluctuation theory of phase transitions is a 
well developed field. In particular, ten thousands of papers in condensed mat- 
ter physics are devoted to the study of pairing fluctuations. At this instant we 
refer to an excellent review of Larkin and Varlamov [15]. 

In some phenomena pairing fluctuations behave similarly to quasiparticles. 
However there are also differences: 



Typical "binding energy" -Ebind of pairing fluctuation is of the same or- 
der as inverse fife-time t^^ oc 1/\T — Tc\, whereas for quasiparticles 



■ Typical size Z oc 1/ ^J\T — Tc\ is large for T near Tc. Quasiparticles 
have a small typical size. 



Fluctuations near behave as classical fields in sense of Rayleigh- 
Jeans: 3 -momentum distribution is n{p) ~ T/E{p) in the vicinity of 



Tc. 



We will demonstrate below that the Ginzburg number {Gi = AT/Tc), 
which determines the broadness of the energy region near the critical tem- 
perature, where fluctuations essentially contribute, is Gi ~ A{Tc/ Hq)'^ with 
A ~ 500 in our case. To compare, for clean metals A ~ 100, fie, the 

latter is the electron chemical potential. Thus Gi 1, if Tc is rather high, 
Tc ~ (| and we expect a broad region of temperatures, where fluctu- 

ation effects might be important. 
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Thermodynamical potential and its mean field solution 

The Fourier component of the density of the thermodynamic potential (ther- 
modynamic potential per unit volume V) in the superconducting quark matter 
with the di-quark pairing can be written in the following form [12, 16], cf. also 
[1,11], 

n = n^ + J2{-co\dr(if + c\V4,\^)+aD + ^D\ (1) 

a,i 

The Greek indices a,/3 = {R,B,G} count colors, the Latin indices i = 
{u, d, s} count flavors. The expansion is presented up to the fourth order in the 
di-quark field operators (related to the gap) assuming the second order phase 
transition, although at zero temperature the transition might be of the first or- 
der, cf. [17]. Clnis the density of the thermodynamic potential of the normal 
state. The order parameter squared is D = jdisP = + |c?Gp + Ids]"^, 
dR\\da\\dB for the isoscalar phase (IS), and D = 3|c/cFLp, Ic^Rp = MgP = 
= IffcFLp, d*j^ ■ da = dc ■ dB = d*^ ■ dR = for the CFL phase, 
da = {d,a,d'^,da}. The so called 2SC phase is the IS phase with unpaired 
quarks of one color (to say R). The interaction of the di-quark field with 
fluctuation gluon fields is usually introduced in the standard way through the 
corresponding gauge-shifted full co-variant derivatives. Recent paper [18] 
demonstrated that the gluon fluctuations contribute essentially to thermody- 
namic quantities only at temperatures in a narrow vicinity of the critical tem- 
perature, AT/Tc ^ 0.1 for relevant values of parameters. In this rather small 
temperature interval, AT, gluon field fluctuations change the nature of the 
phase transition (from the second to the first order). As we will show below, 
the value AT is much less than the temperature region, where fluctuations of 
the di-quark order parameter might be important. Thereby in the latter discus- 
sion we suppress the gluon fields as well as the discussion of any pecuharities 
of this narrow temperature region near Tc. 

Near the critical point coefficients of (1) can be expanded in t = {T — 
Tc) /Tc. In the weak couphng hmit they render: 

a = aoln(T/Tc) ~ aot, ao = ^, (2) 

b, = b, = b/{l+^) = ^-^, C=^^y C0 = SC. 

For the classical (mean) fields 7 = 1 in the IS phase and 7 = 1 /3 in the CFL 
phase, cf. [12], = /ig/S, /x^ is the baryon chemical potential (the contri- 
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bution of the strange quark mass is neglected), C(3) = 1.202... The critical 
temperature Tc is the same for the IS and the CFL phases for 0. For 

nis 7^ 0, Tc would depend on that would result in a smaller Tc for the CFL 
phase than for the IS phase. 

The value of the order^parameter follows by solving the equation of motion 
for the field operators, 5^1 /5d^^ = : 

-cod^dl + cAt4 - a<4 - bDdl = 0. (3) 

The stationary, spatially-homogeneous mean field solution of (3) (without 
taking into account of fluctuations) is 

DMF = -ae{-t)/b, snuF = -^e>{-t), (4) 

where the step function G(-i) = 1 for t < and &{-t) = for t > 0. 

To be specific, discussing T < TcV/e further consider the CFL case. For the 
finite system of a large spherical size R:$> ^, with having the meaning of the 
coherence length, we obtain 



(R 



V2^ 



> ? = WA- (5) 



where A^c = 3 is the number of colors, Nf = 3 is the number of flavors and to 
be specific we assumed the simplest structure j^jp oc 5^. 

Fluctuations of gap in self-consistent Hartree approximation 

Now we will consider fluctuations. Below we present d^^, as = dj^, ^. + 
d^', and above as = d^', where index "c" labels the classical solution. 

Since 5F{V, T) = T), cf. [19], the density of the free energy of the 

CFL phase expanded in (d^')^ terms renders 



SF = dFc + 6F' = J2{-co\drdi^f + c\Vdi^f) + aDc 

k 

k 

+ -H^-^[NcNf (6) 



k 



The last term is introduced within the self-consistent Hartree approximation 
(within the $ functional up to one vertex), uh = 10/9 accounts different 
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coefficients in $ functional for the self-interaction terms (ci^ - for the given 
field d and (4)^(4)^ terms), cf [20]. We presented d^' = Efe <^a,k^''^''''''^ 
k^^ = {u!,k), Xn = {T,—f), and introduced the notation 

E \<k\' = E I'^^i^ = ^ / ^ . '.^ (7) 



(27r)^ cqo;^ — cfc^ — 



In Matsubara technique the temperature dependence is introduced by the re- 
placements: u ^ cOn = 2mnT, -ij^ ^ TEn=-oc/(0- Were- 
stricted ourselves by the self-consistent Hartree approximation. Effects of the 
damping of fluctuations, being produced by the two- and more vertex diagrams 
of are beyond the scope of this simple approximation. As we argue below, 
the main effects we discuss in this paper are not essentially affected by the 
width terms. 

Variation of eq. (5) over yields the spectrum of fluctuations 

cow^ -cP -ni^ = 0, (8) 
with the squared mass parameter 

= r]\a\ = ^{bi + b2/3)D, + a + unih + b2)N,Nf ^ |<^^|2 . (9) 

k 

Notice that the real physical meaning of the effective mass has the quantity 
m/^/cQ rather than m, as follows from (8). 

Solution of the equation for the fluctuating field presented in the coordinate 
space is characterized by the length scale I = by the time scale f = 

^Co^^(c??)"^/2. Above the critical point = 0, = a + 0{J2k \4'k\'^) > 0' 
and neglecting terms one gets the parameter 77 ~ 1. 

Variation of (5) over (i^*^ yields the equation of motion for the classical field 

-(6i + 62/3)L»c<c = 0. (10) 

Only dropping the term responsible for fluctuations we find Dc = Dmf, 
d^a,c = ^L,MF' = "^MF = VmpWI foT T < Tc, cf. (4), (9), and we obtain 
r/MF = 2/3. In reality fluctuations affect the classical solution and renormalize 
the critical temperature of the phase transition Tc yielding TJ*^" < Tc. More 
generally 



m 



2 
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miip + W, 6m^ = --{bi-b2/9)NcNfJ2\'PL\^ (H) 
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for T < Tl^^. From (9), (11) one finds Dc{rn?' = 0) = 0, as the consequence 
of the self-consistency of our approximation scheme. As we have mentioned, 
we neglected fluctuations of the gluon fields, which change the nature of the 
phase transition but yield only a small jump of the order parameter. 

Contribution of fluctuations of gap to specific heat below 

In order to demonstrate how fluctuations of the order parameter may af- 
fect thermodynamic quantities let us calculate the contribution of fluctuations 

to the specific heat density Cy = —T ^^p-j^, cf. [19]. The mean field 
contribution is 

cr = r^e(-i), (12) 

as it follows from (1), (2) and (4). Here and below we use an approximate 

expression a aot valid at T near T^, see (2). The contribution to the specific 
heat from the normal quark excitations is suppressed as oc r~3/2g-A(o)/T 
for T <^ Tc, cf. [21]. We may reproduce correct low temperature behavior 
multiplying (12) by a form-factor / ~ 2e-^^^^/^(Tc/Tf/^il - ^)~^. 

The fluctuation contribution F ' can be determined with the help of the func- 
tional integration 

exp(-(5W) = y^D(/)'exp(-(5W[0']), (13) 
where 5W is an effective potential. Using (13), (5) we obtain 

5F' = -iN.Nf J 0^ln[couj^ - cP - m% (14) 

Again within the Matsubara technique one still should do the replacement 
u^<^n = 27TinT, (0 ^ TE;:=^oo / (gr- We dropped an infi- 
nite constant term in (14). However expression (14) still contains a divergent 
contribution. To remove the regular term that does not depend on the closeness 
to the critical point we find the temperature derivative of (14) (entropy per unit 
volume): 

SS' = ?§^?^N,N,^iKl^. (15, 
The contribution of fluctuations to the specific heat density is then given by 

k 

- -T^N^N.^^g^, ^ -,MPao/T.. (16) 
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In the second line (15) we remained only the terms quadratic in fluctuation 
fields, i.e. we assumed |^^| > |^|. 



Knowing the contribution of fluctuations to the entropy and specific heat we 
may recover their contribution to the free energy and energy densities 

5F' = r S'\p^dT, 5E' = r C{,\p^dT. (17) 
Jo Jo 

With (15) - (17) one may recover fluctuation contribution to all thermodynamic 
quantities. 

Now we are able to calculate the quantity iG {X = 0) = |(/)'^|^ and 
its temperature derivative, where iG (X) is the time-ordered Green function. 
Using the Matsubara replacement lo ^ ojn = 27rinT and the relation J2niy'^ + 
n^)~^ = |cth(7ry), or the corresponding relation between the non-equiUbrium 

Green functions iG {X) and ImG^^* at finite temperature, we arrive at the 
expression 



d^k 1 / \/ ck^ + . 



Using that cth(y/2) = 2nB{y) + 1, where nB{y) = (e^ — 1)~^ are Bose 
occupations, and dropping the regular contribution of quantum fluctuations we 
obtain 

By index "T" we indicate the thermal contribution. 

The integration is performed analytically in the limiting cases. Let m S> 
T^/cq. Then nsiy) — e~^. In the very same approximation one has cP ^ 
m? for typical momenta. Then 

^, ,,2 cl'^ /2mr\3/2 ( m \ ^ 

?i^^i^ = ik(^) ^^^\-T^)^^»^^^^- '''' 

In the opposite limiting case, m ^ TT^fcQ, there are two contributions to the 
integral (19), from the region of typical momenta cfe^ ~ vr? and from the 
region cfc ^ » m^. They can be easily separated, if one calculates the auxihary 

quantity . In the region ck ^ m one may use an approximation 

UBiy) ^ l/y and 

^^M£[cP ~ m'] 0. -(W/2)-^T^. (21) 
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The region ^ m yields a regular term 

12 V c3 ' dT 



k 

General expression for the fluctuation contribution to the specific heat is given 
by the first line (15) and can be resolved with the help of the self-consistent 
solution of (11), (19) (or (20), (21), (22) in the limiting cases). Assuming for 
rough estimates that fluctuations can be described perturbatively and putting 
m? ~ ?Ti^^p = r7MF|a|, from the second line of (15) we find for T <^ 



47r3/2V2cy^|t|V4r2 



and for r> 

87rT2|t|V2 U; c\c) ■ 

In the weak coupling limit for the CFL phase oq/c = 67r2r2/(7C(3)/8), 
Co = 3c, r^MF = 2/3. Eq. (22) holds in the low temperature limit T <^ 
7r(|t|/3)^/2 j"^^ whereas eq. (24) is valid in the opposite limit. The first term in 
(24) dominates over the second one for |t| < 0.7, i.e., in the whole region of 
validity of (24) the main term is the first one. Its singular behavior, ~ 1/ ^/\t\, 
is typical for the specific heat in the vicinity of the critical point of the second 
order phase transition, as in metalhc superconductors. 

Also in the CFL phase (T < Tc) there is a contribution to the specific heat 
of Goldstone-like excitations, cf. phonons in the ordinary condensed matter. 
For T 2> rup^G, where m^.c is the mass of the pseudo-Goldstone excitation, 
we get 

where Nq is the number of pseudo-Goldstone modes. One can see that the 
contribution (25) is numerically small compared to those terms (cf. (22), (24)) 
we have evaluated for temperatures within the fluctuation region. 

Now we are able to discuss the validity of the self-consistent Hartree ap- 
proximation in our problem. Adding the contribution to $ functional with two 
vertices produces the sun-set diagram in the di-quark self-energy. In such a 
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way beyond the Hartree approximation there appears the width term in the 
self-energy (behaving as —i'yoj for small oS). Such a term governs the slow 
relaxation of the order parameter in the phase transition phenomena outside 
the equilibrium. In our case (thermal equilibrium) this term can be dropped 
compared to the cqoo^ term, which we have in the thermodynamic potential (1) 
from the very beginning, at least in both the low and high temperature lim- 
its. Indeed, in the low temperature limit three Green functions entering the 
sun-set diagram produce an extra exponentially small particle occupation fac- 
tor compared to that governs the Hartree term (20). Near Tc, i.e. in the high 
temperature limit, we may remain only n = term in the Matsubara sum 
over frequencies, as we have argued, thus suppressing both the linear and the 
quadratic terms in oj. 

Ginzburg - Levanyuk criterion and Ginzburg number 

Comparing the mean field (12) and the fluctuation (15) contributions to the 
specific heat (in the low and high temperature limiting cases one may use eqs. 
(22), (24)) we may estimate the fluctuation temperature T^^ < Tc, at which 
the contribution of fluctuations of the order parameter becomes to be as impor- 
tant as the mean field one (so called Ginzburg - Levanyuk criterion). 



Fluctuations dominate for T > T^^. For typical values ~ (350-^500) MeV 
and for Tc ^ (50 70) MeV in the weak couphng limit from (26), (22) we 
estimate T^^ ~ (0.6 -i- 0.8)rc. If we took into account the suppression fac- 
tor / of the mean field term oc e~^^^^/'^, a decrease of the mass m due to the 
fluctuation contribution (cf. (11)), and the pseudo-Goldstone contribution (25), 
we would get still smaller value of T^ '^ (< O.STc). We see that fluctuations 
start to contribute at temperatures when one can still use approximate expres- 
sions (22), (20) valid in the low temperature limit. Thus the time (frequency) 
dependence of the fluctuating fields is important in case of CSC. 

In the condensed matter physics one usually performs calculations in the 
high temperature Umit. In this limit one neglects the time (frequency) depen- 
dent terms considering quasi-static thermal fluctuations of the order parameter 
Then the fluctuation contribution is determined with the help of the functional 
integration 




iMF 
V 1 



for T < Tc. 



(26) 



exp(-(5F7r) 

5F'[d'] 




(27) 



i,a,k 
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The integration yields 

5F' = TNfNcJ j^Hck^ + ri\a\]. (28) 

Comparison of eqs. (28) and (14) shows that the high temperature expression 
(28) is recovered, if one drops all the terms except n = in the corresponding 
Matsubara sum over tOn = 27rinT in eq. (14). The contribution of n 7^ 
terms to the J2k l^fcP suppressed in the limit co47r^T^ 3> m^. Then one 
immediately arrives at the specific heat given by the first term of (24). Thus in 
our case one may suppress the frequency dependence of fluctuations only for 
\t\ « 1. 

Simplifying, the energy width of the fluctuation region, where the fluctua- 
tion effects prevail over the mean field ones, is usually estimated following the 
Ginzburg criterion. The probability of the fluctuation in the volume is given 
by ~ exp{-dn{Vii)/T). It is ~ 1 for 6n{Vfi) - T, where Sn{Vii) (the con- 
tribution to the thermodynamic potential in the corresponding variables) is the 
work necessary to prepare the fluctuation within the volume Vq. The minimal 
size of the fluctuation region characterized by an order parameter d ^/Dmf 
is ^(r)vW»7' cf. (5), (8). Thus, taking Sn{Va) = 5r2MF(Vii) = for the 
typical temperature T^, when fluctuations start to dominate, we obtain 

^ - 2b 3 - 3 ■ ^ ^ 

Although the above estimate is very rough we took care of all the numeri- 
cal factors. This allows us to notice that the value Tq^^ estimated from (26), 

if one uses eq. (24) for Cy remaining there only the first term, is ~ 
0.5|i(r|;<)/t(rgG)|i/2rG for = Nf = 3. All the dependencies on the 
parameters in expressions for T^^ and were proven to be essentially the 
same. 

Fluctuation region is determined by the Ginzburg number Gi = \Tc — 

Tq\/Tc. The larger the value Gi is the broader is the fluctuation region. For 
clean conventional superconductors [15] Gi ~ A{Tc/ fi)^ ^ 10~^^ 10~^^, 
A ~ 80, whereas for superfluid He^ and in our case Gi ^ 1 (since A ~ 500 
and T/iJ.q can be as large as ^ | in favorable cases). The role of fluctuations 
increases in cases, when the effective dimensionality of the droplet decreases 
or/and when the quark mean free path decreases. Both possibilities result in 
a significant increase of the Gz-number. In the case of quark droplets of the 
typical size L <C ^ one deals with the zero-dimensional system, Gi oc $,^/L^ 
and the fluctuation contribution to the specific heat behaves as Cy ~ for 
Gi <C \t\ <C 1. For dirty superconductors Gi increases oc l/{pFeie)^ where 
PFe is the electron Fermi momentum and Ig is the electron mean free path. A 
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significant decrease of the quark mean free path is expected due to the presence 
of the hadron impurities inside the quark droplet. Thus, there is still a variety 
of possibilities for a further enhancement of superconducting fluctuations. 

Fluctuations of the gap above 

All above expressions for fluctuating quantities (except the vanishing of the 
pseudo-Goldstone contribution above Tc and the appearance of an extra di- 
quark decay contribution to the width) are also valid for T > Tc, if one puts 
Dc = in general expressions or m^^p = a > 0, = 0, rj = t^mf = 1 
in the corresponding approximate expressions. Above Tc we should compare 
the fluctuation contribution to a thermodynamic quantity with the quark and 
gluon contributions of the normal state of the quark-gluon plasma. In the weak 
coupling limit for the specific heat density of the quark-gluon plasma one has, 
cf. [22], 

The first two terms are quark contributions and the third term is the gluon 
contribution. We omitted a contribution of sttange quarks which is rather small 
for T < rus and we neglected the temperature dependent effective gluon mass. 
We also disregarded the corrections to the quark and gluon terms since 
such corrections were not taken into account for condensate quantities. If we 
included the quark-gluon masses that appear in the framework of the approach 
[23] matching the quasiparticle quark-gluon description and the lattice results, 
we would get even smaller contribution of Cy, that is in favor of fluctuation 
effects. 

Above Tc, eqs. (22) and (24) yield rather smooth functions of T except 
the region t <^ 1. The fluctuation region is rather wide since (5y ~ (5y 
even at T essentially larger then T^. The appearance of an extra channel of 
the di-quark decay width, 27^60'^, beyond the Hartree approximation does not 
qualitatively change the situation since 7dec is a regular function of T, 7dec on 
T — Tc, and vanishes in the critical point. Thus the applicability of the high 
temperature limit (one can put aj„ = in the Matsubara sum) is preserved in a 
wide temperature region, T ~ Tc, [15]. Only far above Tc situation might be 
changed. The main uncertainty comes from the values of coefficients (2) which 
were derived for T rather near Tc and in the weak coupling limit. Bearing all 
this in mind we estimate the value of the fluctuation temperature T^'^ ~ 2Tc 
for typical value Tc///^ ~ 0.1 ^ 0.3. The higher this ratio is, the stronger is the 
contribution of fluctuations at the given ratio T/Tc. Note that a decrease of fiq 
with increase of the temperature results in an increase of the ratio Td iXq{T). 
On the other hand jjLq should be at least larger than zero temperature gap in 
order one could speak of any pairing fluctuations at given temperature.. 
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Thus we see ihai fluctuations of the di-quark gap may essentially contribute 
to the thermodynamic quantities even well below Tc and above Tc. 

Assumptions which we have done 

Note that several simplifying assumptions were done. The coefficients (2) 
were derived for |i| <C 1 in the weak coupling limit neglecting the strange 
quark mass, but apphed in a wider temperature region for may be not suffi- 
ciently large jj^g. Only Gaussian fluctuations were taken into account within 
the self-consistent Hartree approximation. Thus, di-quark width effects were 
assumed to be suppressed. We used the expansion of the thermodynamic po- 
tential up to the fourth order terms in the mean field neglecting a small jump 
in the order parameter due to gluon fluctuations. As was argued in [18] gluon 
fluctuations may also contribute but in a more narrow vicinity of Tc than the 
order parameter fluctuations (for the values of the ratio Td iiq ^ 0.15 0.3, 
which we are interested here). We incorporated in the thermodynamic quan- 
tities only the terms which have a tendency to an irregular behavior near Tc, 
whereas above Tc the short range correlations begin to be more and more im- 
portant with the increase of the temperature. Therefore one certainly should 
be cautious applying above rough estimates outside the region of their quan- 
titative validity. However, as we know from the experience of the condensate 
matter physics, see [15], such extrapolation equations work usually not too bad 
even for temperatures well below and above Tc. 

Fluctuations of temperature, density, magnetic susceptibility 

So far we have discussed the specific behavior of fluctuations of the order 
parameter at fixed temperature and density. There are also fluctuations of the 
temperature and the local quark density. They are statistically independent 
quantities [19], < ST6pq >= 0, and their mean squares are 

The averaging is done over the volume, P is the pressure, Va is as above the 
volume related to the fluctuation. In the limiting case Cy <^ Cy, we ob- 
tain < {STf > /T^ ~ pq/{N^C^^), where A^^ is the number of quarks 
involved in the volume Vq_. Thus, far from the critical point the contribution 
of fluctuations is suppressed as \/< (ST)'^ > jT'^ oc 1/ N^- This standard 

fluctuation behavior is essentially changed for Cy ^ Cy. Then we get even 
larger suppression of the temperature fluctuations. For |t| <C 1 in the high 
temperature limit we obtain ^< (ST)'^ > /T^ ^ \t\ — > for |t| — > 0, where 
we assumed that fluctuations are most probable within the typical fluctuation 
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volume Vfl = 47r^^/3 and ^ oc l/vl^I for t ^ 0. Thus, if the system is at the 
temperature T in a narrow vicinity of Tc, fluctuations of the temperature are 
significantly suppressed. Being formed at < — t ^ 1 the condensate region 
evolves very slowly, since the heat transport is then delayed (the typical evolu- 
tion time is roughly r oc 1/ At temperatures T outside a narrow vicinity 
of Tc, fluctuations of the temperature resulting in a significant decrease of the 
temperature in the volume Vq are rather probable, < (5T)^ > /T^ ~ 1. 
This is the consequence of a very short coherence length and, thus, not too 
large number of quarks contained in the volume Vg, ^ 2± 0.13/(rcy^) being 
~ 0.5 fm for Tc ~ 50 MeV and for |t| ~ 1. Fluctuations of the quark density 
are also large for typical \t\ ^ 1, due to the smallness of Thereby, we argue 
that the system may produce the di-quark condensate regions of the typical size 
^, thus feeling the possibility of the phase transition even if its temperature and 
density are in average rather far from the critical values. 

Other quantities associated with second derivatives of the thermodynamic 
potential are also enhanced near the critical point demonstrating typical 1 / 
behavior, cf. [21]. However numerical coefficients depend strongly on what 
quantity is studied. E.g. fluctuation contributions above Tc to the color dia- 
magnetic susceptibilities 

Xa = -(d^SF/dHl) Ha = cm\Aa, (32) 

are proven to be ^ 1 everywhere except very narrow vicinity of the critical 
point. In spite of a smallness, as we know, in metals the fluctuation diamag- 
netism turns out to be of the order of the Pauh paramagnetism even far from 
the transition. Also for T S> Tc contribution of fluctuations to the magneto- 
conductivity of 2D electron systems is experimentally distinguishable, [15]. 

How gap fluctuations may manifest in heavy ion collisions 

Anomalous behavior of fluctuations might manifest itself in the event-by- 
event analysis of the heavy ion collision data. In small (L) size systems, L < ^, 
(zero dimension case would be L ^ the contribution of fluctuations of the 
order parameter to the specific heat is still increased, as we have mentioned, 
see [15]. The anomalous behavior of the specific heat may affect the heat 
transport. Also kinetic coefficients are substantially affected by fluctuations 
due to the shortening of the particle mean free paths, as the consequence of 
their rescatterings on di-quark fluctuations. If thermalization happened at an 
initial heavy ion collision stage, and a large size system expands rather slowly, 
its evolution is governed by the approximately constant value of the entropy. 
Due to a large contribution to the specific heat (and to the entropy) of di-quark 
fluctuations, an extra decrease of the temperature may occur resulting in an 
essential slowing of the fireball expansion process (due to smaller pressure). 
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Having the di-quark quantum numbers, fluctuations of the gap may affect the 
di-lepton production rate from the quark-gluon plasma, as it has been noted in 
[13]. Another question is how one can distinguish di-quark fluctuations related 
to the CSC from the quark fluctuations, which may relate to the deconfinement 
transition? Theoretically, since is, in general, different from the deconfine- 
ment transition temperature T^ec, in assumption that Tc > T^ec there might 
be two anomalous fluctuation peaks related to two distinct collision energies. 
E.g., if the observed peak [25] in the i^+/7r+ ratio at ~ (30 40) GeV/A is, 
indeed, associated with the deconfinement transition [26], it would be worth to 
seach a possibility of another peak at a somewhat higher collision energy than 
the first one, related to the CSC phase transition, and vice versa, if the observed 
peak relates to the CSC, another peak at a somewhat lower energy could relate 
to the deconfinement. However it remains unclear is it possible to distinguish 
these peaks experimentally, if Tc and T^ec are rather close to each other. One 
needs a careful measurement of /7r~^ ratio in the whole collision energy 
interval. 

Pairing fluctuations in hybrid stars 

Besides the crust and the hadron shell, the hybrid star contains also a quark 
core. Both the nucleon shell and the quark core can be in superconduct- 
ing phases, in dependence on the value of the temperature. Fluctuations af- 
fect transport coefficients, specific heat, emissivity, masses of low-lying ex- 
citations and respectively electromagnetic properties of the star, Hke electro- 
conductivity and magnetic field structure, e.g., renormalizing critical values 
of the magnetic field (Hd, He, Hc2)- Note that the CSC is the type 1 super- 
conductor with the Ginzburg-Landau parameter k, < 0.6, if < 14 MeV 
(estimation is done for //g ~ 400 MeV), and it is the type 11 superconductor 
with K essentially larger than one, if Tc can indeed be large, Tc <i 50 MeV, that 
is the case which we are interested in here, cf. [18]. Especially the type IT CSC 
could have interesting consequences for observations, cf. [11]. 

The effect of thermal pion fluctuations on the specific heat and the neu- 
trino emissivity of neutron stars was discussed in [27, 28] together with other 
in-medium effects, see also reviews [29, 30]. Neutron pair breaking and for- 
mation (PBF) neutrino process on the neutral current was studied in [31, 32] 
for the hadron matter. Also ref. [32] added the proton PBF process in the 
hadron matter and correlation processes, and ref. [33] included quark PBF 
processes in quark matter. PBF processes were studied by two different meth- 
ods; with the help of Bogolubov transformation for the fermion wave function 
[31, 33] and within Schwinger-Kadanoff-Baym-Keldysh formalism for non- 
equihbrium normal and anomalous fermion Green functions [32, 28, 29]. 
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As was observed in [32, 28], analogously to the Direct Urea (DU) process 
n — ^ pei>, these processes (n nvv and p pvv) have a large (one-nucleon) 
phase space volume, if the pairing gap is A > 1 MeV in some density interval 
(in addition to that DU and PBF emissivities have similar exponential suppres- 
sion factors oc e~^/-^). Moreover, these processes are affected by nucleon- 
nucleon [32] and electron-electron [34, 35, 30] correlations in such a way 
that the emissivity of the process n nvv is not significantly changed but 
the emissivity of the process p pfv increases up to ten - hundred times. 
This enhancement is due to the fact that the square of the bare vertex for the 
p pvu reaction contains a very small ~ ~ 0.006 factor compared to 
the corresponding factor ~ 1 for the n —>■ nvv reaction channel. However 
the proton may produce the neutron-neutron hole by the strong interaction and 
the electron-electron hole by the electromagnetic interaction, which then may 
couple to the weak current. This circumstance is still ignored in a number of 
works, which rediscovered this process with vacuum vertices and used it within 
the cooling code, e.g., see [36]. It seems that an artificial a ten - hundred times 
suppression of the relevant process may affect their conclusions. Numerical 
simulation of the neutron star coohng that incorporated PBF processes with 
inclusion of correlation effects, as well as other relevant in-medium effects 
of the nucleon-nucleon and nucleon-pion interaction, like softening of pion 
modes, was performed in [37]. The PBF processes in the quark matter are also 
affected by correlation effects but, as for the case of the reaction n nvv, in 
the given case the correlation effects are expected to be not so significant. 

Order parameter fluctuations allow for extra neutrino processes having no 
exponential suppression in a broad region of temperatures near Tc. 

Contribution of pairing fluctuations to the specific heat in the hadron shell 
is minor for the case of the neutron pairing due to a small value of Tc ^ IMeV 
compared to the value of the neutron chemical potential > 50 MeV). 
Therefore in the neutron channel fluctuations of the gap are relevant only in a 
very narrow vicinity of the critical point. However this effect might be not so 
small for protons, for which the chemical potential is of the order of several 
MeV, whereas the gap is of the order of one MeV. Therefore it seems that 
fluctuations may smear the phase transition in a rather broad vicinity of the 
critical point of the proton superconductivity. 

Pairing fluctuation effects in the quark matter are especially important, if the 
critical temperature of the CSC phase transition is rather high {Tc ^ 50 MeV), 
as is estimated for 2SC and CFL cases. Temperatures ~ 10 ^ 50 MeV may 
indeed arise at an initial stage of the hybrid star cooling. The fluctuation con- 
tribution to the total specific heat for T < Tc is evaluated in a line with that we 
have done above. As noted in ref. [15], the Aslomasov - Larkin contribution 
to the electro-conductivity proves to be the dominant term in the vicinity of 
Tc. Analogously, one may expect a significant increase of the heat conductiv- 
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ity and viscosity that would result in a delay of the heat transport at the initial 
stage of the hybrid star cooling. Moreover, neutrino may efficiently rescatter 
on fluctuations that governs their drift to the hadron shell at T < T^. These 
effects may be especially important for the CFL phase, since in that case all 
other mechanisms of the heat transport are suppressed at T < Tc. Estimation 
of all these effects needs however a separate study. 

Concluding remarks 

We would like once more to emphasize that the coefficients of the thermo- 
dynamic potential (1) obtained in the weak coupling limit might be essentially 
modified, in case if we applied the results to the fireball produced in a heavy 
ion colUsion or to hybrid stars, somehow changing our conclusions. E.g., the 
quark chemical potential decreases with the temperature. It results in an ad- 
ditional increase of the contribution of fluctuations, cf. dependences of the 
coefficients (2) on jiq. Thus we may need to analyze the strong coupling limit 
instead of the weak coupling limit we discussed. There exist arguments that 
the strange quark mass is very large and due to that the phase transition 
in the CFL phase does not occur up to very high baryon densities [38]. If Hq 
becomes smaller than the strange quark mass we come from the possible 3SC 
phases to the 2SC phases. A more general discussion of modifications beyond 
the framework of the weak coupling limit can be found in [24]. A discussion 
how the coefficients (1) may in general vary within the Ginzburg - Landau 
approach is given in [12]. Our above arguments are admittedly speculative, 
and mean only to demonstrate a qualitative possibility of their application to 
heavy ion collisions and to an initial stage evolution of the hybrid star. A more 
detailed treatment of the problem definitely needs a further work. 

Concluding, within the Ginzburg - Landau approach we estimated the fluc- 
tuation energy region at the CSC phase transition. The quantitative estimates 
are based on the values of parameters derived in the weak coupling limit for 
fluctuations of the CFL order parameter. Qualitative results survive also for 
fluctuations of other possible phases. We found that the frequency dependence 
of fluctuations is important for CSC in a wide temperature region. Fluctua- 
tions may contribute essentially to the specific heat even at T rather far be- 
low and above Tc. We estimated Tq^< < O.STc and rfl,> ~ 2Tc. Our rough 
estimates show that the high temperature CSC could manifest itself through 
fluctuations of the di-quark gap in the course of the heavy ion collisions at 
SIS 300, if the critical temperature of the phase transition is indeed rather high, 
Tc ^ (50 70) MeV. The CSC fluctuations are also relevant for an initial stage 
of the hybrid star evolution, if Tc is ~ 50 MeV. However quantitative results 
depend on the values of several not sufficiently known parameters and further 
studies are still needed to arrive at definite conclusions. 
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